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A procedure for using the eigenvectors of the elements of the representations of a discrete group in 
model building is introduced and is used to construct a model that produces a large reactor mixing 
angle, sin 2 813 = | sin 2 jg, in agreement with recent neutrino oscillation observations. The model 
fully constrains the neutrino mass ratios and predicts normal hierarchy with the light neutrino mass, 
mi « 25 meV. Motivated by the model, a new mixing ansatz is postulated which predicts all the 
mixing angles within la errors. 



Introduction: We use the group SU(3) and its discrete 
subgroup S*4 for model building. Both of these groups 
had been studied extensively as flavour symmetry groups, 
e.g. [1] [2]. The S4 group has the presentation 

(a,b\a 2 = b 3 = (ab) 4 = e). (1) 

54 is the symmetry group of the cube (Fig. 1) and the 
elements of the group can be represented as the orien- 
tation preserving rotations of the cube. The matrices 
representing the generators can be written as 

/ooi\ , /non . . 

0=0-10, 0=100. 2 
Vi 0/ Vo 1 0/ w 

Here the basis vectors e\ — (1,0, 0) T , e-i = (0,1, 0) T 
and e 3 = (0, 0, 1) T form the symmetry axes of the cube 
passing through face centres. If we define the left-handed 
leptons, L = (L e ,L fl ,L T ) T where L e = (v e L,eL) T etc., 
as a triplet in this basis, the flavour states L e , L^, L T 
correspond to e\, e^, e% respectively. Usually in models a 
set of flavons are introduced whose vacuum expectation 
values (vevs) produce the desired texture for the fermion 
mass matrices. In other words, the orientation of fermion 
flavour states as well as the flavon vevs in the flavour 
space determines the form of the mass matrices. 
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FIG. 1. The generators o and b represent 7r-rotation about 
axis a and ^-rotation about axist, respectively. 

Axes of the orientation preserving rotations of the cube 
are nothing but eigenvectors of the corresponding rota- 
tion matrices with eigenvalue equal to +1. The basis 



vectors e%, e2 and e 3 are examples. There are also other 
vectors like the ones passing through the opposite edge 
centres (e.g. axis a in Fig. 1) and the ones passing through 
the opposite vertices (e.g. axisf, in Fig. 1). Compared to 
vectors pointing in random directions, these vectors are 
"special" in the context of the S4 symmetry. The ro- 
tation matrices are unitary and so their eigenvalues in 
general are complex numbers with unit modulus. If non- 
degenerate, these eigenvalues also correspond to unique 
eigenvectors and the author argues that these eigenvec- 
tors are also "special" like the rotation axes. 

As an example consider v = (l,w,w) T , the normalised 
(v'v — 1) eigenvector of the matrix b in Eq. (2), cor- 
responding to the eigenvalue u), where u = e 1 ^ and 
ui = e~ l ~£~ . Since e l0 v, where e 10 is an arbitrary phase, 
is also a normalised eigenvector, we impose the following 
condition to uniquely fix the phase: The component of 
the eigenvector in the direction of one of the basis vectors 
should have zero phase ie. &rg(v' ej) = 0, where i = 1, 2 or 
3. This is intuitive since the basis vectors are used to de- 
fine the fermion flavour states. Thus the allowed choices 
for v are ^=(1, ui,w) T , ^g(w, l,w) T and 4=(cD,a;, 1) T . 

Let g represent an element of the group and I be one of 
its non-degenerate eigenvalues. The corresponding "spe- 
cial" eigenvector, eig(g, is defined using the following 
normalisation condition and the phase condition: 

eig(fl, /)] eig(s, i)i = 1, arg (eig(a, l)\ e^J = 0. (3) 

For example, the basis vector is eig(a&, 1)3. Such eigen- 
vectors which will be used later in the model are listed 
below: 



eig(6, l)i = 




eig 


;0, 1)1 = 


73 


eig(6,w)x = 




ei 


g(c,i)i = 


^(i,o, ? ) T , 


eig(6,w)i = 


^j(l,w,w) T , 


ei 


g(c,i)3 = 


^Ho,i) T 


eig(d, l)i = 




eip 


(c,-«)i = 


£(1,0, -if 


eig(d, 1)3 = 


^(-i,o,i) T , 


eip 


(c,-«) 3 = 





(4) 

where c = bob and d = aibab) 2 . 



The Model: Recent experiments [3] have shown that the 
reactor mixing angle, #13, is non-zero. A number of mod- 
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els and parameterisations have been proposed, e.g. [4], to 
accommodate the non-zero #13. The model described in 
this paper is constructed in the Standard Model (SM) 
framework with the addition of the right-handed neu- 
trino triplet, v R — (yiR,V2R,V3R) T in the context of the 
type-1 seesaw mechanism. We postulate a global flavour 
symmetry group, 



the symmetry breaking, the charged-lepton mass term, 
Eq. (7), takes the form 



l L T^M d l R + H.C 



where 



1 1 1 



Gf = SU(3)i x SU(3) 2 x U{l) f . 



(5) 



v3 \i u q. 



(9) 



(10) 



The fermion fields and a set of postulated flavons belong 
to specific representations of G/ as shown in Table I. 
The U(l) group is introduced to ensure that the flavons 
couple to only the desired fermions. We write the mass 
terms at the lowest order, containing the fermions and 
the minimum number of flavons, invariant under G / and 
the SM gauge group. The eigenvectors of the elements 
of the S*4 subgroup of the SU{3) group are used to con- 
struct the vevs of the flavons. The vevs break the flavour 
symmetry and the required mass matrices are obtained. 



(e R ,fj, R ,T R y and M d = 
= etc. The charged- 



h = (e L ,/i L ,r L ) T , l R - 
diag(m e , m M , m T ) with m e 

lepton mass matrix, T' M<i, when left- multiplied with T, 
is diagonalised giving the charged-lepton masses m e , m^, 
m T . T is the Trimaximal mixing matrix [6] . For later use 
we define 2x2 maximal matrices, B2 and B3: 



B-2 = 





(11) 
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U(l) f 


fl 


fl + fe 


fi + U 


fl+fr 







-u 







fl 






For the neutrinos, the lowest order Dirac mass term is 



H.C. 



(12) 



TABLE I. The fields <j> e , (j)^, <f> T , (f>, ^1, £2 are the flavons. 
For the 1/(1) group the tabulated values are the generators, 
e.g. ft = e lf ' e . The SM Higgs is a flavour singlet. 



The tensor product expansion of the fundamental rep- 
resentations of SU(3) are 



where H is the conjugate Higgs and y w is the coupling. 
The flavon £1 belongs to the representation 6 of SU(S)i 
and can be written as a 3 x 3 complex-symmetric ma- 
trix. The restriction of 6 of SU(3) to the S4 subgroup 
is the direct sum of 1, 2 and 3 of S4. We assign a very 
simple choice of vev for £1 where 2 and 3 parts vanish, 
ie. (£1) becomes the identity when written in the ma- 
trix form. After the symmetry breaking, the Dirac mass 
term, Eq. (12), takes the form 



6. 



6. 



3(83 = 108. (6) 



m w [Wvr + {vrY{v l ) c \ +H.C 



(13) 



For the charged leptons, the lowest order mass term is where Ul — [v t v l v l) 1 



L (y e (f>ee R + y M <t>^R + vAtT R ) H-^- + H.C. (7) 

where y a are the coupling constants, H is the SM Higgs, 
A c j >l is the cut-off scale for the flavons </> e , (f>a and (j> T . 
The ^4 group has four irreducible representations: 1, 2, 
3 and 3' [5]. The orientation preserving rotations of the 
cube discussed earlier belong to 3'. The flavon triplets 
4> e , 4>n an d 4>t belong to the representation 3 of SU(3)i. 
The restriction of the representation 3 of SU(3) to its 
subgroup S4 is the representation 3' of S4. Therefore 
the "special" eigenvectors of the representation matrices 
of 3' of 5*4 are used to construct the vevs of the flavons. 
We assign: 

(<t>e) = eig(M)i, (<M = ci gO> w )i: (<M = eig(6,w)i 

where the angular brackets are used to denote vevs. The 
flavon vevs break the flavour symmetry and the Higgs 
vev, (0,h o ) T , breaks the weak gauge symmetry. After 



The lowest order Majorana mass term for the neutrinos 



is 



yG{vR) C ^24> VR.J—J? 



H.C. 



(14) 



Note that <fi transforms as a 3 under both SU(3)i and 
SU(3)2- Therefore (j) can be written as a 3 x 3 matrix, 
<j>ij, the row index i representing SU(3)i and the column 
index j representing SU(3)2- The flavon £ 2 belongs to 
the representation 6 of SU (3)2- A 6 of SU(3), just like a 
6, contains a 1, a 2 and a 3 of the S4 subgroup. As was 
done earlier for the case of (£1), here we assign (£ 2 ) also 
to be equal to the identity. After the symmetry breaking, 
the Majorana mass term, Eq. (14), takes the form 

m G (^F(0) {4>) T vr + H.C. (15) 
where vtiq = . VG . a . The matrix Uf))Uf)) T is complex- 

symmetric and (cf>)((f>) T contains all the interesting 
physics in our model. 
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To assign vev for the flavon cfnj, we use (S^i x (84)2, 
the subgroup of SU (3)i x SU(3) 2 . The group (5*4)1 x 
(5*4)2 nas 24 x 24 = 576 elements. Let gi and g 2 be the 
elements of (5*4)1 and (54)2 respectively. If v\ and v 2 
are the eigenvectors of 171 and 172 corresponding to the 
eigenvalues a\ and a 2 , then the direct product v\ x v 2 
will be an eigenvector of g\ x g 2 with an eigenvalue (Zict2. 
Now we make the following assumption: 

((/)) = v\xv\ + vl X vl + v\xvl + vfxv% (16) 

where the RHS of Eq. (16) is the sum of four eigenvectors. 
Based on the choices for v\ x v 2 s we get a set of similar 
cases of solutions described in the following sections. The 
assumed form of (</>) given in Eq. (16) and the choices for 
v\ x ujs were obtained through educated guesses and also 
through trial and error to fit the experimental data. 
Case 1: Here we assign 



(17) 



"l 1 


= e 2 


> 


vl 


= e 2 




= cig 


'-.(a, l)i, 


vl 


= ei 




= ci§ 


;(<M)i, 


vl 


= ei 


vi 


= eig 


id,i) 3 , 


vl 


= ei{ 



Using Eq. (16), Eqs. (17) and Eqs. (4), we get 

\ + \/2 



{<t>) = 



(18) 



in matrix form, where the row and the column indices of 
the matrix correspond to the (54)1 and the (54)2 indices 
respectively. If ma » m w , the Majorana mass matrix, 
rriG{4>){<t>) T from Eq. (15), becomes much larger than the 
Dirac mass matrix, m w x Identity from Eq. (13), resulting 
in the type-1 see-saw mechanism. We get an effective see- 
saw mass matrix [7], M SSl proportional to the inverse of 

<0><0) T : 



-k((4,)(<t>y 



where k — 



(19) 



M. 



(20) 



is diagonalised using the unitary matrix U v 
Uj,M ss U* = -diag (mi, m 2 ,m 3 ) 
where the neutrino masses mi, m 2 , m 3 are given by 

fc(2 + V2) fc(2 + V2) 



m 1 



and 



with 



m 2 



m 3 



1 + J2{2 + y/2) 



2(2 + V2) 
(21) 



U v = B 2 XB 2 EB T 2 V 

X = diag(l,l,i), 

£ = diag(e l f,l,l), 

V = diag(e~ 4 ra,l,e^ra). 



(22) 

(23) 
(24) 
(25) 



The PMNS matrix becomes 

U = TU V = TB 2 1B 2 £B^V. 



(26) 



The mixing obtained is a constrained form of the 
Trichimaximal (T%M) mixing [8] with x = yg : 



1^1 = |TXM (X _ 



(27) 



The modulus sign is used throughout this paper to indi- 
cate that the expression for the mixing matrix is valid 
only upto right and left multiplication with diagonal 
phase matrices (which do not affect the phenomenon 
of neutrino oscillation). The right multiplying diagonal 
phase matrices, like V in Eq. (26), do contribute to Ma- 
jorana phases, the study of which is beyond the scope of 
this paper. From Eq. (27) and using Eq. (12) in [8] we 
get 



|t/ e3 | 2 = " sin 2 ^ =► sin 2 013 = 0.025, 
16 



\U e - 



sin 2 0i2 = 0.342, 



sin 2 023 



S C p = 



1 

2' 



(28) 
(29) 
(30) 
(31) 



From Eq. (21) we get the ratios of the neutrino masses, 
mi : m 2 : m 3 = 0.945 : 1 : 2.117. These ratios are 
compatible with the mass-squared differences measured 
experimentally [9, 10] within la errors and thus we can 
predict the light neutrino mass: 



24.7 meV < m x < 25.5 meV. 



(32) 



Even though the set of eigenvectors given in Eqs. (17) 
results in the matrix (<p) given in Eq. (18), other choices 
of eigenvectors also exist which produce the same (<j>). 
The matrix T t (^)(^) T T* is a highly constrained form of 
the complex-symmetric "Simplest" texture [11]. 

Let x = i n where n is an integer and let 



$ = 

1 T 




(33) 



Using Eq. (18) can be rewritten as 

(<(>)= K (34) 

Case 2: Assigning v\ = e 2 , v\ — e 2 , v\ = eig(c, i)i, 
vl = ei, vf = eig(c, z) 3 , t| = ei, vf = eig(d, l)i and 
v\ = eig(c,-i) 3 we get 

(0)=$* 2 . (35) 
In this case, the resulting PMNS matrix is 

U = TU V = TB 2 1 2 B 2 £B^V (36) 
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and it is a constrained form of the Triphimaximal (T0M) 
mixing [8] with (f> — — jg : 

|CT| = |T0M ( ^_ ft) |. (37) 
From Eq. (37) and using Eq. (9) in [8] we get 

2 ' ' / 2 T ~) sin 2 023-0.613, (38) 



where 



l^3| 2 -gSin 2 



2tt 



SCP = 7T- 



(39) 



Here, as well as in the next two cases, Eq. (21) which 
gives the neutrino masses and Eqs. (28, 29) which give 
sin 2 #i3 and sin 2 812 remain valid. 

Case 3: Assigning v\ — e%, v\ — v\ — eig(a, l)i, 



eig(d, 1)3, «| = ei, vf = eig(d,l)i and 



v\ — eig(c, -i) 1 we get (<ft) — $* 3 resulting in 

U = TB 2 1 3 B 2 £B^V => \U\ = |T X M (x= _ „ } |, (40) 

' (41) 



sin 2 823 



2 



TCP 



3tt 
Y' 



(42) 



Case 4: Assigning v\ = e 2 , v\ — e 2 , v\ = eig(c,-z)i, 
eig(c,-i)3, w| = ex, vf = eig(d, 1) 3 and 



ex, v\ = 



v\ = eig(c, -i) 3 we get (<f>) = <E>* 4 resulting in 



U = TB 2 2 4 B 2 £B 2 r V 



sm 



2tt tt 
T + 16 



|C/| = |T0M (0=ft) |, (43) 

=> sin 2 6*23 = 0.387, (44) 
: 2tt. (45) 



The values predicted by all the four cases of the 
model are within 3cr errors of the experimental best fits 
[9, 10, 12]. In fact the generic prediction sin 2 6x3 — 0.025, 
Eq. (28), is within la errors. However the global analy- 
sis [9] shows more than 2a tension with sin 2 823 = \ , the 
T%M value (Cases 1 and 3, Eqs. (30, 41)). On the other 
hand the T^M values, sin 2 8 23 = 0.613 from Eq. (38) in 
Case 2 and sin 2 823 = 0.387 from Eq. (44) in Case 4, are 
well within la errors calculated in [10] and [9] respec- 
tively. All the cases predict sin 2 8 12 = 0.342, Eq. (29), 
which is at the edge of the 2a error range in [9]. A new 
mixing ansatz called the VS mixing [13] is proposed in 
the following section which modifies 8x2 as well as 5qp- 
The VS Mixing Ansatz: The mixing obtained using 
the model, Eqs. (26, 36, 40, 43), is of the form 



\U\ = \TB 2 T l B 2 £Bl 



(46) 



The matrix \TB2l n \ gives the Tribimaximal (TBM) mix- 
ing [14]. Multiplying TB 2 l n with B 2 £Bj mixes the first 
and the third columns of the TBM matrix. Similarly we 
may also mix the first and the second columns resulting 
in the new ansatz defined by 



|VSi»(a)| - \TB 2 T l B 2 £BlB3E'Bl 



(47) 



£' = diag(e M ,l,l). 



(48) 



Note that the Cases 1 to 4 are simply VSjn (0) with n 
1 to 4 respectively. Eq. (47) on simplification gives 



where 



|VSin(a)| = \TB 2 T l H2SH 3 \ 



5 = diag(e i T6,l,l) 










is\ 






is' 




n 2 = 





1 


I 




IS 


d 






\is 
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with 



a 



16' 



16' 



(49) 

(50) 
(51) 

(52) 



We get sin 2 6x2 within la errors for 0.087T < a < 0.267T. 
Table 2 lists a few cases of the VS mixing along with the 
predicted values of the mixing angles. The author finds 
the choice a = | to be aesthetically pleasing. When 
a = I we get £' = £ and also d — c, s' = s. 





VS.i(f) 


VSi(f) 


VS.i(f) 


VSi(f) 


VS-i(f) 


VSi(f) 


sin 2 623 


0.613 


0.387 


0.613 


0.387 


0.613 


0.387 


sin 2 612 


0.323 


0.323 


0.317 


0.317 


0.319 


0.319 


8 c p 


1.277T 


0.27tt 


1.187T 


0.187T 


1.137T 


0.137T 



TABLE II. Note that sin 2 613 = 0.025 is a generic feature of 
the VS mixing. Conjugation, VS*n(a), changes the sign of 
5c p without affecting the mixing angles 612, &23 and 6*13. 

Summary: The symmetries represented by a discrete 
group are related to the eigenvectors of the group ele- 
ments. We develop a notation to uniquely identify the 
eigenvectors and use it to assign vevs for the flavons. 
An orthonormal set of eigenvectors define the fermions' 
flavour states. The model thus constructed predicts the 
reactor mixing angle, sin 2 8x3 — 0.025, as well as the ra- 
tios of the neutrino masses, mi : m 2 : m 3 = 0.945 : 1 : 
2.117. The T(f>M versions of the model provide solutions 
for 823 in the first octant, sin 2 6*23 = 0.387, as well as in 
the second octant, sin 2 823 = 0.613. The T0M as well as 
the T%M versions give sin 2 8x2 — 0.342. A new mixing 
ansatz, VS,n(a), is introduced which gives reduced val- 
ues for 8x2- The ansatz also predicts various values for 
Scp- 
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